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Extending the scheme developed for a single mode of the electromagnetic field in the preceding 
paper "Structure of multiphoton quantum optics. I. Canonical formalism and homo- 
dyne squeezed states", we introduce two-mode nonlinear canonical transformations depending 
on two heterodyne mixing angles. They are defined in terms of hermitian nonlinear functions that 
realize heterodyne superpositions of conjugate quadratures of bipartite systems. The canonical 
transformations diagonalize a class of Hamiltonians describing non degenerate and degenerate mul- 
tiphoton processes. We determine the coherent states associated to the canonical transformations, 
which generalize the non degenerate two-photon squeezed states. Such heterodyne multiphoton 
squeezed are defined as the simultaneous eigenstates of the transformed, coupled annihilation op- 
erators. They are generated by nonlinear unitary evolutions acting on two-mode squeezed states. 
They are non Gaussian, highly non classical, entangled states. For a quadratic nonlinearity the 
heterodyne multiphoton squeezed states define two-mode cubic phase states. The statistical prop- 
erties of these states can be widely adjusted by tuning the heterodyne mixing angles, the phases 
of the nonlinear couplings, as well as the strength of the nonlinearity. For quadratic nonlinearity, 
we study the higher-order contributions to the susceptibility in nonlinear media and we suggest 
possible experimental realizations of multiphoton conversion processes generating the cubic-phase 
heterodyne squeezed states. 

PACS numbers: 03.65.-w, 42.50.-p, 42.50.Dv, 42.50.Ar 



I. INTRODUCTION 

In this paper we extend to bipartite systems of pairs 
of correlated modes the multiphoton canonical formalism 
developed for a single mode of the electromagnetic field 
in the companion paper "Structure of multiphoton 
quantum optics. I. Canonical formalism and ho- 
modyne squeezed states" which, from now on, we 
will refer to as Part I. Preliminary to the description of 
the multiphoton canonical formalism we need to discuss 
briefly the well established formalism of two-mode quan- 
tum optics 12, 13 • Two-mode non degenerate squeezed 
states are generated by second order susceptibility con- 
tributions excited by laser shots on nonlinear optical me- 
dia. Due to the second order nonlinearity, the frequency 
f2 of a high energy laser splits inside the crystal into a 
pair of frequencies uji,uj2 {O. = wi -\- UJ2) associated to 
correlated modes of the electromagnetic field (non de- 
generate down conversion process), with the reversed up 
conversion process of recomposition of the two frequen- 
cies being allowed as well. The degenerate limit loi = uj2 
can be considered, with the splitted frequencies being as- 
sociated to a single mode of the electromagnetic field. 
Two-photon squeezed states are the coherent states of 
the Hamiltonian which describes the simultaneous cre- 
ation or annihilation of a photon in each of the corre- 



* Electronic address: 'dellanno@sa.in tn.itl 
t Electronic address: desiena@sa.infn.it] 
■f Electronic address: >illumi nati@sa.in fn.itl 



lated modes (or of two photon in a single mode for the 
degenerate case). The two-photon Hamiltonian is di- 
agonalized by coupled linear canonical transformations. 
The two-photon squeezed states can be obtained as com- 
mon eigenstates of the two-mode transformed variables. 
They can be generated as well by applying a (squeez- 
ing) unitary operator on the two-mode vacuum or on the 
two-mode coherent state; the unitary operator moves the 
original mode operators to the canonically transformed 
variables. 

When considering possible extensions of the above 
scheme to multiphoton processes, one natural choice is to 
consider down-conversion processes in which the carrying 
laser frequency splits into n frequencies, giving rise to the 
creation (and later annihilation) of n photons in n corre- 
lated modes of the field (or in a single mode of the field in 
the degenerate limit). However, the related Hamiltonians 
cannot be exactly diagonalized, and an exact canonical 
scheme cannot be implemented. In this paper we follow 
a different route, outlined in Part I for single-mode sys- 
tems. We consider classes of multiphoton processes in 
multimode systems such that the number of photons in- 
volved in the processes does not correspond in general to 
the number of excited modes. In such instances, the fully 
degenerate limit entails that the interaction terms in the 
Hamiltonians do not in general reduce to simple pow- 
ers of the creation and annihilation operators. In Part I 
0, we have introduced degenerate, canonical n-photon 
generalizations of the degenerate, canonical two-photon 
processes. The single-mode multiphoton generalization 
has been obtained by using nonlinear operator functions 
of homodyne superpositions of conjugate quadratures. 
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To the end of generalizing the canonical formalism 
to two-mode, non degenerate multiphoton processes, we 
consider nonlinear operator functions of heterodyne su- 
perpositions of the two-mode variables, obtaining cou- 
pled two-mode transformed operators with a nonlinear 
dependence on the original operators. As in the one- 
mode instance, we show that the canonical conditions 
reduce to simple algebraic constraints only on the com- 
plex coefficients of the transformations, while they are 
independent of the form of the nonlinear functions. 

Upon introducing multiphoton canonical transforma- 
tions for bipartite systems, we analyze the problems re- 
garding the multiphoton Hamiltonians induced by the 
transformations, their physical interest, and their ex- 
perimental realizability. We determine the multiphoton 
coherent states associated to the transformations, and 
study their statistical properties. Finally, we show that 
the two-mode multiphoton coherent states can be gener- 
ated by acting with a unitary operator on the two-mode 
vacuum. The structure of the two-mode multiphoton 
canonical transformations does not allow the realization 
of the associated coherent states by the displacement op- 
erator method, i.e. by the unitary evolution with the 
interaction Hamiltonian. This problem is similar to that 
of mapping interacting fields onto free fields in Quantum 
Field Theory for nonlinear systems. We thus consider 
the coherent states defined as the simultaneous eigen- 
states of the coupled canonical transformations. Exploit- 
ing the entangled state representation we can then com- 
pute their eigenfunctions. The resulting states share im- 
portant features of the standard coherent states, such as 
forming an overcomplete basis in Hilbert space. These 
multiphoton coherent states associated to the nonlinear 
transformations share squeezing properties as well, and 
we thus name them (two-mode) heterodyne multipho- 
ton squeezed states (HEMPSS). They are non Gaussian, 
highly non classical, entangled states. A large part of 
the present paper is devoted to the characterization and 
quantification of their physical properties. 

The independence of the canonical transformations by 
the form of the nonlinear functions allows in principle to 
introduce an infinite number of diagonalizable Hamilto- 
nians. In analogy to the degenerate case we have con- 
sidered the subclass of Hamiltonians involving only alge- 
braic powers of the heterodyne variables. This choice de- 
fines Hamiltonians associated with degenerate and degen- 
erate multiphoton processes in two modes of the electro- 
magnetic field. The higher-order multiphoton processes 
are related to higher-order contributions to the suscepti- 
bility in nonlinear media. We then study in detail only 
the simplest case of quadratic nonlinearity, which leads 
to Hamiltonians describing up to four-photon processes. 
We analyze the photon statistics of the HEMPSS, show- 
ing the strong dependence of the interference in phase 
space by the competing phases entering in the nonlinear- 
ity, and by the heterodyne angles. We suggest schemes 
for a possible experimental realization of the HEMPSS 
by exploiting contributions to the susceptibility up to 



the fifth order in nonlinear media. We determine the 
unitary operator acting on the vacuum to generate the 
HEMPSS: its form suggests alternative routes to their 
possible experimental realization. 

We remark that multimode nonlinear canonical trans- 
formations of the form introduced in the present paper 
can be defined in principle for generic bosonic systems 
of atomic, molecular, and condensed matter physics. In 
such cases however the constraints on the standard form 
of the kinetic energy term and on the conservation of the 
total number of particles play a crucial role and must be 
taken carefully into account. A thorough analysis of the 
relevance of the nonlinear canonical formalisms for ma- 
terial systems lies thus outside the scope of the present 
paper. 

The paper is organized as follows. In Section ^] we 
introduce the two-mode nonlinear canonical transforma- 
tions and the associated multiphoton Hamiltonians. In 
Sections IIIII and Hvl we review the entangled state repre- 
sentation, and exploit it to determine the coherent states 
of the transformations, defined as the eigenstates of the 
transformed annihilation operators (heterodyne multi- 
photon squeezed states, or HEMPSS). In Section W\ we 
derive the form of the unitary operators generating the 
HEMPSS by acting on the two-mode vacuum. In Sec- 
tion we study the photon number distribution of the 
two-mode multiphoton squeezed states. In Section IVIll 
we review the theory of quantized fields in nonlinear me- 
dia, and in Section rVIIII we present a proposal for a possi- 
ble experimental generation of the HEMPSS. Finally, in 
Section IIXI we draw our conclusions and discuss further 
future developments. 



II. TWO-MODE MULTIPHOTON CANONICAL 
TRANSFORMATIONS VIA HETERODYNE 
PHOTOCURRENT VARIABLES 

In this Section we introduce the two-mode multiphoton 
canonical formalism. We consider two (correlated) modes 
of the quantized electromagnetic field. The modes are 
characterized by the annihilation and creation operators 
oi, 02, a\ and at^, obeying the canonical commutation 
relations 

[ai,a]] = 6ij , [ai,a.j] = . (1) 

We begin by recalling that the form of the two-mode, 
linear canonical squeezing transformations is Q: 

bi — ^ai + va^ , 

(2) 

62 = /ia2 -I- vaX , 

with V complex numbers. As is well known, the above 
transformations are canonical if 

= (3) 
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and the constraint is automatically satisfied by the 
paramctrization 

fi = coshr , v = sinhre"^ , (4) 

where r is the squeezing parameter. The transformations 
Eqs. 101 can be obtained from the original mode opera- 
tors ai, 02 by acting with the unitary operator Q 

5i2(C) = exp(C4aI - C*«2ai) , C = re'^ • (5) 

The transformed variables bi , 62 define the diagonalized 
Hamiltonian 

H = b\bi + h\b2 , (6) 

which can be written in terms of the fundamental mode 
operators as: 

H = + W\^)ia\ai + ala2) + 2\iy\^ 

+ 2^*va\a\ + 2^v*aia2 ■ (7) 

The scheme defined by relations Eqs. ijU, 10), de- 
scribes non degenerate two-photon down-conversion pro- 
cesses. 

We wish to extend the linear canonical scheme to 
provide an Hamiltonian description of multiphoton pro- 
cesses. Let us define a new set of transformed variables 
hi and &2 by 

(8) 

62 = ^a2 + va\+xF\a\,a2) ■ 

Here /i, z/, 7 and x s-re complex numbers. The generaliza- 
tion of the two-photon linear canonical transformations 
is obtained by introducing a nonlinear, operator-valued 
function F of the creation and annihilation operators. 
The function F must satisfy some regularity conditions. 
Even so, the canonical conditions may not in general be 
satisfied by a completely arbitrary, though well behaved, 
function. Wc will now introduce two crucial characteri- 
zations of F, allowing to select a large class of functions 
that do satisfy the canonical constraints. First, we as- 
sume the argument of F to be of the form 

aa\ + I3*a2 , a, P e C . (9) 

We next require that 

F^{Z)=F{Z^). (10) 

The bosonic canonical commutation relations that must 
be satisfied by the transformed modes are 

h,b]]=S,, [6„6,]=0. (11) 

We now show that they are satisfied if some simple al- 
gebraic constraints hold between the coefficients of the 



nonlinear transformations. These algebraic relations gen- 
eralize Eq. (PJ of the linear case, and are independent of 
the particular form of the function F. The first condition 
is again Eq. J^Jl. It derives from the c- number part of the 
canonical commutation relations Eqs. (|11() which, in turn, 
is due to the linear part of the transformations Eqs. JHJ. 
The new conditions derive from the operator-valued part 
of Eqs. Hll|l . Exploiting the constraints on the nonlinear 
function F given by Eqs. (|5))- (fTn)) . and using the commu- 
tation rules for functions of the annihilation and creation 
operators (see also Part 10), one finally has 

\a\'-\(3f=Q, 

(12) 

MX*" - i^X*/3* + M*7/3* - = . 

We see that conditions Eqs. (|12|l are simple algebraic 
relations on the coefficients of the transformations. Con- 
ditions Eqs. l|^- H12|) are compatible with very many spe- 
cific canonical realizations of the nonlinear transforma- 
tions. To restrict further the possible canonical choices, 
we will adopt the parameterizations Eq. Q and 

7 = \l\e''' , X=|x|e^'^ , 

(13) 

a = |a|e*^i , p = |/3|e**^ . 

Let us finally choose \a\ = = so that the moduli 
(strengths) I7I and |x| of the nonlinear couplings, and 
the phases Si, S2, 0i, and 62 are the remaining free pa- 
rameters. The argument of the nonlinear function F in 
Eq. © becomes (e"'^=a2 + e'^^a\)/V2, and it can be 
interpreted as the output photocurrent of an ideal het- 
erodyne detector 0, . We see that the structure of the 
two-mode, multiphoton canonical transformations is in- 
trinsically based on heterodyne mixings of the different 
field quadratures. The transformations Eqs. ((SJ now read 

e_j.2^j^\ 

(14) 

52^.a. + .at + |,|e-^-^( --""^^-""^ ), 

while Eqs. ^ are encoded in the single condition 

coshrlxle^'^'^^"^^) - sinhr|x|e"*(''^+^^"'^^ 
-hcoshr|7|e*(''i"''^)-sinhr|7|e''(*i+''i-'^)=0. 

(15) 

Eq. H15|l can be written in the compact form 
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tanhr 



(lxP + l7ncos(0i+i 



+ 2|x||7l cos(5i + (52 - 0) + - |7P) M^i + ^2 - 



IxP + |7P + 2|x||7l cos((5i + (52 + 01 + ^?2 - 2(/.) 

I 



(16) 



Being r real, the imaginary part of Eq. (|16|l must be set 
to zero: 



(17) 



(|xp-|7nsin(0i+02-0) = O. 
Eq. H17|) is satisfied either if 

9i+92-(l)^kTr, (18) 

or 

1x1 = l7|. (19) 
By choosing Eq. H18() then Eq. (|16|l reiiuces to 

tanhr = ±1, (20) 



which correspon(is to r 
Eq. l(Ti1)l we have 



oo. By choosing instea(i 



tanh r 



3s(6>i +02 -<!>) + cos((5i + 62-4 
1 + cos((5i + i52 + 6'i + 6*2 - 20) 



(21) 



which is independent of the nonUnear strength I7I and 
can always be satisfied. The class of possible solutions 
of Eq. H21|l is rather large. In order to find explicit ex- 
amples we must specify some concrete realizations. We 
then adopt the limiting solutions: 



^51 



0,±7r , 01 



±7r,0. 



(22) 



The transformations Eqs. (|14|l define a diagonalized 
Hamiltonian of the form Eq. JB)). Inserting in this ex- 
pression Eqs. (|14|l for the transformed variables yields 
multiphoton Hamiltonians written in terms of the funda- 
mental variables a.i and al. They are parametrized by F, 
and each one is characterized by a specific choice of the 
nonlinear function. Among all the many possible forms 
we select some which bear particularly interesting and 
realistic physical interpretations. Let us in fact consider 
Hamiltonians associated to the following set of nonlinear 
functions defined as F{() = C", i.e. as integer powers 
of the heterodyne variable. They describe nondegener- 
ate and degenerate processes up to 2n-photon ones. In 
particular, we concentrate our attention on the simplest 
choice of lowest nonlinearity F{C) = C^, describing up to 
four-photon processes. The Hamiltonian reads 



H = Aq + _Bo(a|ai + 0^02) 
+Co{a\^ + (^2 '^2 + 2a|aia2a2) 



-[-Dia|a| + D2a'ia2 + -D2'^ia2 

,t2„ J , „t„t2 



t2 



-|-L'5(a{ axo"^ -f ix^a'2 0.2) + h.c. 



(23) 



where the coefficients are given by: 

Ao = I7P + 2IH', 

Bo = UA' + W\' + 2\j\\ 

Co 



1 



D' = 



1 



|7|e2'«i(e"'>* +e~*"V), 



iS2 , 



^3 = ^|7|e'''^(e'^>*+e"^^V), 



h\ 21(81+62) 



_ I ,|2„i(ei-i-e2) 



, ^5 = l7re 



• (24) 



The parameters are constrained by the canonical condi- 
tions Eqs. H22|) . Regarding Ao, Bo,Co, it is sufficient to 
insert the parametrization Eqs. Q). To be concrete, for 
the remaining parameters we take the canonical choice 
(5i -I- (52 — = and 9i + 62 — (j) = it, obtaining 

Di =e^'^(sinh2r-2|7p). 



D2 = 



hi^r^i{2ei~52) ^/ ^ _hl^r^i{2e2-5i) 
2 2 



D4 = 



|7P 



2id 



2 id 



(25) 



As previously stated, the Hamiltonian Eq. (|23|l describes 
up to four-photon degenerate and non degenerate pro- 
cesses. We conclude this Section noting that for degen- 
erate, single-mode processes, the heterodyne canonical 
formalism reduces to the homodyne canonical formalism 
introduced in Part I 0. 



III. ENTANGLED STATE REPRESENTATION 

In this and in the following Sections we determine the 
form and the statistical properties of the coherent states 
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associated to the heterodyne canonical transformations, 
at least in the simplest cases. In order to do this, it is con- 
venient to introduce preliminarily the so called entanqled 
v? "T? '77 Oi 

state representation 6, 7, 8, tfj. 

We define the non Hermitian operators : 



Z 



0,2 



P. 



V2 ^/2 
which satisfy the commutation relations 



(26) 



(27) 



Note that Z is just the argument of the nonlinear func- 
tion F in the transformations Eqs. H14|l : as already men- 
tioned, this operator can be interpreted as the output 
photocurrent of an ideal heterodyne detector 0, |^ . 
In terms of Z, P^, Eqs. (|14|l become 



bi 



z^ + i 



fi' - V" 



P, + |7|e^*'iF(Z), 



(28) 



hi - 



%/2 



Z + 1- 



V2 



Pi + \iW'-F{Z^) , 



where ^l' = e*^i/i, ^J^" = e'^^^, v' = e^^'^ii', v" = e"'''^;^. 
The algebra defined by relations Eqs. (|27l) and 12t)|l is 
characterized by the orthonormal eigenvectors {\z)} of 
Xq^ + Xg^ and — Pg^ [H , also called entangled-state 
representation, where Xg. and Pg. = Xg^+iL, i = 1,2, are 
the homodyne quadrature operators for the mode a,, and 
z is an arbitrary complex number z = zi + iz2. It can be 
easily proved that in the two-mode Fock space the state 
\z) is 



exp 



exp 



+ V22 



\piz* a 



al at 



100) 



\/2za], 



V2, 



z ag^ 



exp 



-at at 



|00),(29) 



where |00) is the two-mode vacuum. The states Eq. H29f) 
satisfy the eigenvalue equations 

{Xg^+Xg^)\z) = 2zx\z) ,{Pg^~Pg,)\z) ^2Z2\Z) , (30) 

and 

Z\z) ^ z\z) Z^\z)^z*\z). (31) 
The states {\z)} are orthonormal 

(z'|z) =^,5(2)(/-z), (32) 
and satisfy the completeness relation 

J d'^z\z){z\ = 
- /d^z : exp{-2|zp + V2z(at + agj 

TT J 2 

+V2z*{al^ + aej - (aj^ + agj{al^ -t- aej} : = 1 , 

(33) 



where (Pz = dzdz*. In the next Section we will exploit 
the entangled state representation to determine the co- 
herent states associated to the nonlinear canonical trans- 
formations. 



IV. HEMPSS: HETERODYNE MULTIPHOTON 
SQUEEZED STATES 

As is well known, the Glauber (harmonic oscillator) 
coherent states can be defined via three equivalent pro- 
cedures ■ For Hamiltonians whose elements belong to 
more complex algebras than the Weyl-Heisenberg alge- 
bra, the three procedures are in general not equivalent, 
and lead to different definitions of coherent states 10, 111. 
The coherent states associated to the unitary Hamilto- 
nian evolution applied to the ground state are generated 
by acting with the displacement operator on the vacuum. 
This is not our case: as anticipated in the Introduction, 
we have been faced with the question of constructing op- 
erator variables which form a harmonic-oscillator Weyl- 
Heiscnberg algebra, although containing nonlinear terms. 
Having solved the problem in the entangled state rep- 
resentation, we can now determine the coherent states 
defined as simultaneous eigenvectors of the transformed 
operators Eqs. (|14|) . The associated eigenvalue equations 
are 

= /3,|V')/3 , i = l,2, (34) 

which in the {\z)} representation become 

{z\b,\i^)0 = . (35) 

Exploiting Eqs. 128|) and the differential representation 
for the operators {Z,Pz}, we can write Eqs. (|35|l in the 
differential form 

^' ^ ""-zi^p + ^^-^d,,i>[3 + |7|e^^ip(z*)V/3 = /3i# : 



V2 



V2 



(36) 



%/2 V2 
whose solution is 

ll:p{z, Z*) = Ne-<^\A''+Tiz'+T,z-B{z,.') (37^ 

In Eq. H37() . iV is a normalization factor, the function 
B[z, z*) (which encodes the nonlinear contribution) is 

B{z, z*) = 6i y" dCFiC) + b2 f d^FiO , (38) 



and 



/i' + u" _ + v' _l + 2iIm{^jL'*v") 

iS2 



/i' - v" /i" - v' 



bi = 



16 1 



\/2|7|e 
/i'-i/" ' 

V2/3i 



bo - 



%/2|7|e 
/i" - I/' ' 

V2p2 



(39) 
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We remark that the canonical conditions Eqs. (|22|l imply 
Re[a] > 0, /r7i[a] — and Re[B{z, z*)] — 0, which ensure 
the normalizability of the wave function Eq. (|37|l . The 
states Eq. (|37|l satisfy the (over)completeness relation 



d^Pid^f]2Wp{i;\ = 1. 



(40) 



Of course, for 7 = the states Eq. H37|l reduce to the 
standard two-mode squeezed states. The states Eq. (|57|) 
are two- mode multiphoton states, depending on hetero- 
dyne combinations of the field modes. They share prop- 
erties of coherence and squeezing, as we will show in 
more detail in the following. We thus name them het- 
erodyne multiphoton squeezed states (HEMPSS). For a 
single mode, they reduce to the homodyne multiphoton 
squeezed states (HOMPSS) introduced in the compan- 
ion paper Part I. We have explicitly obtained the general 
form of the wave functions of the HEMPSS in the entan- 
gled state representation. The general properties of these 
states, which we will study in the following Sections, are 
however easier to investigate in the coordinate represen- 
tation. _Let_ us first recall the following useful expression 
for \z) 



_ g — 2iziZ2 



(41) 



where the tensor product (g) \(j>) denotes kets in the 
Hilbert space Hi (E) H2, and is an eigenvector of 

the quadrature operator Xg. of the z-th mode {i — 1,2). 
From Eq. l|lT|l it follows 



where tpfj is evaluated at zi — ^^^^^"^ . The states 
Eq. (|42ll are in rather implicit form. Under particular 
conditions one can however find explicit analytic expres- 
sions. For instance, considering a quadratic nonlinearity 
F{Z) = Z"^, and letting Si — 61 = ^, they take the form 



N 



■ exp{- 



[{xe^ - 3:9 J + (Fi - Fz; 
4a — 6iS(2;ej -I- xg^) 



r „fXg-^ -f Xg^ 3 i^Bi + Xg^ 2 

• exp{-z^( Y - a( ) 

+(rr + F.)(:^^)}, 



(43) 



V. UNITARY OPERATORS 

The heterodyne multiphoton squeezed states 
(HEMPSS) defined in the preceding Section are 
not the states unitarily evolved by the Hamiltonians 
associated to the nonlinear canonical transformations. 
As in the degenerate case, however, they can be unitarily 
generated from the two-mode vacuum. By looking at 
the entangled state representation of the HEMPSS, 
Eq. l|H7|) . it is evident that 



|V')/3 = U{Z., Z^)Di{ai)D2{a2)Si2{g)m 



(44) 



where 



^CK^a 1 — a ■ ai 



(45) 



with g — ±r and ai = ^i* Pi — v[3* (i ^ j — 1,2). These 
two operators generate a two-mode squeezed state from 
the two-mode vacuum. The operator U takes the form 

U{Z, Zt) = exp I -61 j diF{i) -hj d^F(e)| , 

... 

and the canonical conditions assure that it is unitary. 
Obviously, for 7 = 0, the state |V')/3 reduces to an or- 
dinary two-mode squeezed state. In the case of lowest 
nonlinearity F{Z) ~ Z^, exploiting the canonical condi- 
tions, the operator Eq. H4t)|) reads 



U{Z, Z^) = e 



(47) 



dz2e2*("''2~"«i)"^V/^(2i,Z2), (42) where 



V2| 



7|e 



3(/i 



In this expression the exponent of the operator Eq. H47|) is 
clearly anti-hermitian. Depending on the specific canoni- 
cal choice, the parameter r can acquire both positive and 
negative values. It is interesting to express the operator 
Eq. I|47|) in terms of the rotated modes ag. and a\:. 



U = exp 



A 

2V2 



«ei + SaJ^ae^ + SaJ^Og^ + a^J + h.c. 



,t2 



(49) 

In Section rVIIII wc will use this form to propose a possible 
experimental realization of the HEMPSS. As it is evident 
by their expressions, the HEMPSS are non Gaussian, en- 
tangled states of bipartite systems. 



where, due to the canonical conditions, 

2%/2, , coshr - sinhre*(^i+* 
-I7I- 



3 cosh 2r — sinh 2r cos(0i -I- 62 
is a real number. 



VI. PHOTON STATISTICS 

In this Section we study the statistics of the four- 
photon HEMPSS. We first compute the two-mode pho- 
ton number distribution (PND). We compare it with the 



7 



PND of the two-mode squeezed states, thoroughly stud- 
ied in Ref. i and we show its strong dependence on the 
modulus and phases of the nonlinear couplings, and on 
the heterodyne mixing angles. We next show the behav- 
ior of the average photon numbers in the two modes as 
functions of |7|. We then look at the second order corre- 
lation functions, whose behaviors are strongly dependent 
on the parameters of the nonlinear transformations as 
well. 

Let us consider the scalar product {ni,n2\z) 



.^™i-m^n2-m^(Af-m)(2|zn . 



771 ! 

M! 



1/2 



(50) 



Here m = 777177.(711,77.2), M = 77iax(7ii, 712) and L"(x) are 
the generalized Laguerre polynomials. Eq. 15U|I . together 
with the relation of completeness Eq. 1331) , allows to write 
the following integral expression for the photon number 
distribution (PND): 



P(71l,7l2) = 1(711, 712 I lA)/?!^ = 



7n! 
M! 



(51) 



d^ze-l^l z*"i-"z"=^-™i(^'^-'")(2|znV'^(z) 



This expression can be computed numerically, and plot- 
ted for different values of the parameters. We compare 
the four-photon HEMPSS with the standard two-mode 
squeezed states. The case 7 = yields the PND com- 
puted in Ref. • In Fig. ^ the plot is drawn for an 
intermediate value of the squeezing parameter (r = 0.8), 
and a low mean number of photons, while in Fig. [3 we 
take a strong value of the squeezing parameter (r = 1.5) 
at a larger mean number of photons. The PND is sym- 
metric with respect to iii and 712, and exhibits the char- 
acteristic oscillations due to interference in phase space 
|l3l |: the oscillations are enhanced for higher squeezing. 
We have plotted the two-mode, four-photon PND for 
different values of the parameters. In all graphs we have 
kept the same values of r, and the corresponding val- 
ues of Pi , (32 used for the standard squeezed states 
(r = 0.8 , /5i = /32 = 3 ; r = 1,5 , /?i = /32 = 5). 
We first set to zero the heterodyne mixing angles and 
the phase of the squeezing: 6*1 = 6*2 = = 0. This 
choice forces Si + S2 = tt (see Eq. (HU). Among all the 
possible realizations of this constraint we have selected 
the completely symmetric choice 61 — S2 = 7r/2, and the 
completely asymmetric one 61 — n , 62 — 0- For each 
value of I7I we have then four plots for different values of 
the squeezing and of the phases. We choose, respectively, 
I7I = 0.1 and I7I = 0.2. 

We see that the form and the oscillations of the PND 
for the four-photon HEMPSS strongly depend on the 
values of the phases Si and S2 which produce compet- 
ing effects in U{Z, ^t). In fact, in Figs. El [7| El cor- 
responding to the symmetric choice for the phases, the 




FIG. 1: P{ni,n2) for a two-mode squeezed state, for r = 0.8, 
/3i = /32 = 3 and I7I = 0. 




FIG. 2: P(ni,n2) for two-mode squeezed state, for r — 1.5, 
/3i = /32 = 5 and I7I = 0. 



PND is symmetric with respect to 711 and 712. Viceversa, 
in Figs. El EllHl the asymmetric choice of the phases 
leads to an asymmetric PND: the peaks are displaced, 
and their number and form are changed. 



As for the standard squeezed states, enhancing the 
squeezing parameter r enhances the oscillations of the 
PND. The effect of the nonlinear strength I7I is compet- 
ing with that of r. In fact, we see that in Eq. H48|l the 
effective strength A of the nonlinear term involves the 
product |7|e~'", with the choice of parameters yielding 
r > 0. Then, a sufhciently high value of r reduces the 
effect of increasing |7|. If we look now at the symmetric 
PND's (Figs. El El El EI), we see that, for a high value of 
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FIG. 3: P(ni,n2) of the four-photon HEMPSS, for r = 0.8, 

/3i = /32 = 3 and I7I = 0.1, 91=02^ 0, Si ^ S2 = f . 




FIG. 4: P(ni,n2) of the four-photon HEMPSS, for r = 0.8, 
/3i = /32 = 3 and I7I = 0.1, 61=62^ 0, (Si = tt, (Sa = 0. 



r, increasing I7I from to 0.2 does not change apprecia- 
bly the oscillations (see Figs.|Sl|51coinpared with Fig.jSJ. 
On the contrary, if r is fixed to a lower value, increasing 
I7I results in a smoothening of the oscillations (Figs. 13 
[3 compared with Fig. Obviously, when I7I attains 
very high values, it cannot be contrasted by r, and the 
oscillations are even more suppressed. In Figs.llllandll2l 
we have plotted the PND for /3i = P2 = 0, and for differ- 
ent values of r; this corresponds to the two-mode four- 
photon squeezed vacuum. As in the case of the standard 
two-mode squeezed vacuum we obtain a diagonal, 
symmetric PND, in spite of the unbalanced choice on the 
phases. 

In order to study the influence of the local oscillators 
angles, in Fig. ^| we have plotted the PND for the same 
parameters of unbalanced case of Fig. |S1 but with di = 




FIG. 5: P{ni,n2) of the four-photon HEMPSS, for r = 1.5, 
/3i = /32 = 5 and I7I = 0.1, 61=82 = 0, 81=62 = f . 




FIG. 6: P{ni,n2) of the four-photon HEMPSS, for r = 1.5, 

Pi=P2 = 5 and I7I = 0.1, 01=92 = 0, <5i =n,52 = 0. 



—62 = 7r/6; we see that the PND is modified with respect 
to that of Fig. |S1 with a slight re-balancing between the 
two axes, and less pronounced peaks. 

In Figs. I14landll5lwe show the mean numbers of pho- 
tons in the two modes as a function of I7I for different 
values of the other parameters. We see that the balanced 
choice on the phases gives < ni >=< 712 > (correspond- 
ing to a symmetric PND), and the mean photon num- 
bers increase monotonically with increasing |7|. In the 
unbalanced instance (corresponding to a non symmet- 
ric PND), the average photon number in the first mode 
is markedly larger and monotonically increasing, while 
the average photon number in the second mode first de- 
creases, and increases monotonically beyond a certain 
value of I7I. This behavior is similar to that exhibited 
by the HOMPSS for single-mode systems, as discussed 
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FIG. 7: P{ni,n2) of the four-photon HEMPSS, for r = 0.8, 
/3i = /32 = 3 and |7| = 0.2, 61=62= 0, 5i = 82 = f . 




FIG. 8: P{ni,n2) of the four-photon HEMPSS, for r = 0.8, 
/3i = /32 = 3 and I7I = 0.2, 61=62= 0, <5i =11,62= 0. 

in Part I Q. 

We next study the average photon number for the 
HEMPSS as a function of the local oscillators angles 9i, 
02, comparing with the case of two-mode squeezed states 
fFig. ll6|l . To this aim, when I7I = we retain as a formal 
trick the dependence on 9i and 62 by fixing (f) — 9i + 62, 
according to canonical conditions Eqs. 1221) . Because of 
the choice /3i = /32, the plot (ni) = (712) is somehow 
redundant; it is symmetric and shows an oscillatory be- 
havior, as it should be. 

For the nonlinear case, we fix the parameters in the 
following way: we let /3i = /32 = 3, r = 0.8 , I7I — 0.4 
and Si — T^, 62 = 4>, because of the canonical conditions. 
As we can see in Figs. El and the situation is very 
different from that of the standard two-mode squeezed 
states. The shape is strongly deformed and we can ob- 
serve regions with a suppressed or enhanced number of 




FIG. 9: P{ni,n2) of the four-photon HEMPSS, for r = 1.5, 
f3i=f32 = 5 and {j] = 0.2, 61=62 = 0, <5i = <52 = f . 




FIG. 10: P(ni,n2) of the four-photon HEMPSS, for r = 1.5, 
Pi=P2 = 5 and I7I = 0.2, 61=62 = 0, <5i = n, 82 = 0. 



photons with respect to the mean number of a standard 
two- mode squeezed state (see Fig. II6II . 

Finally, we have plotted the second-order two-mode 
correlation function (/'^■'(O) = (o|aia2a2)/(a|ai) (0202) 
as a function of the local oscillators angles 61 and 02- 
We again obtain a variety of shapes, corresponding to 
different statistical behaviors. In particular we observe 
oscillating transitions from bunching to antibunching be- 
haviors. In this case as well, we show for comparison 
the three-dimensional plot of ^(^^(O) for a two-mode 
squeezed state Fig. ^1 It is of course symmetric and 
exhibits a narrow region of antibunching. 

Already for small strengths of the nonlinearity (I7I = 
0.1) we can see from Fig.EOlthat the shape of the correla- 
tion function is deformed with respect to that of the two- 
mode squeezed state. This deformation becomes more 
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FIG. 11; P(ni, 712) of the four-photon HEMPSS, for r = 0.8, 2){ni) = (712) for 61=62 = ^ (dotted line). 
I3i= 132 = and I7I = 0.1, 61=62 = 0, Si =71,82 = 0. 




HEMPSS, at r = 1.5, /3i = /32 = 5, Si = 02 = 0: l)(ni) 
(full line) and (712) (dot-dashed line) for (5i = n, S2 = 0; 



FIG. 12: P(ni,n2) of the four-photon HEMPSS, for r = 1.5, 2){ni) = (712) for (5i = 52 = f (dotted line). 
/3i = /32 = and |7[ = 0.2, 61 = 62 = 0, 5i =71,^2= 0. 




FIG. 13: P(ni,n2) of the four-photon HEMPSS, for r = 0.8, 
I3i=l32 = 3 and I7I = 0.2, 61 = -62 = f , <5i = tt, ^2 = 0. 



FIG. 16: (ni) = (712) as a function of 61 and 62 (0 = +^2), 
of a standard two-mode squeezed state (7 = 0), for r = 0.8, 
01=02= 3. 
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FIG. 17: (ni) as a function of 61 and O2, of the four-photon 
HEMPSS, for r = 0.8, /3i = /32 = 3, I7I = 0.4, 5i = vr, (Ja = 0. 




FIG. 18: (na) as a function of 0i and 62, for the four-photon 
HEMPSS, at r = 0.8, /3i = /Ja = 3, I7I = 0.4, 5i ^11,62 = 4>. 




FIG. 19: ff^^>(0) as a function of 9^ and 6*2(0 = 6*1 + 6^2), for 
a standard two-mode squeezed state with 7 = 0, r = 0.8, 
Pi =132 = 3. 




FIG. 20: ff<2'(0) 

as a function of 6*1 and 02, for the four- photon 
HEMPSS, at r = 0.8, Pi = P2 = 3, |7i = 0.1, h = tt, and 
52 = (p- 



(HI <0) 




FIG. 21: ff'^HO) 

as a function of 61 and 62, for the four- photon 
HEMPSS, at r = 1.5, Pi = P2 = 5, I7I = 0.1, Si = tt, and 
52 = 4>. 



marked for larger values of I7I, as shown in Figs. 1221 and 
From Fig. 12] we see that for increasing values of the 
squeezing parameter, the antibunching behavior is atten- 
uated and the effect of the nonlinearity is progressively 
suppressed. This is a further indication of the competing 
effects of the squeezing and of the nonlinearity. 

Summing up, similarly the case of single-mode sys- 
tems, we observe that the photon statistics of the two- 
mode, four-photon HEMPSS can be strongly modified 
by the gauging of the parameters entering the nonlinear 
canonical transformations, such as the heterodyne mixing 
angles 9i and 62, the strength I7I, and the phases Si and 
62 of the nonlinearity. Let us notice that this flexibility 
can be further enhanced by considering the large number 
of different possible canonical choices for the parameters 
that have not been explored in the present work. 
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FIG. 22: g^^^O) as a function of 6*1 and 62, for the four-photon 
HEMPSS, at r = 0.8, (3i = P2 = 3, I7I = 0.2, 5i = tt, and 
62 = 4>. 




FIG. 23: ^(^'(O) 

as a function of 6*1 and 62, for the four-photon 
HEMPSS, with r = 0.8, /3i = /Ja = 3, lij = 0.2, Si = tt, and 

($2 = 



VII. QUANTUM DYNAMICS IN NONLINEAR 
MEDIA AND GENERATION OF THE HEMPSS 

The interaction part of the Hamihonian Eq. H23|) con- 
tains many non degenerate and degenerate multiphoton 
terms; this complex structure is evidently due to the non 
trivial canonical constraints associated to the structure of 
the nonlinear Bogoliubov transformations. On the other 
hand, the very existence of a canonical structure suggests 
that the experimental realization of the associated mul- 
tiphoton processes and squeezed states should be con- 
ceivable. To this aim, we must bear in mind that any 
experimental realization of the HEMPSS must necessar- 
ily involve high order contributions to the susceptibility 
in nonlinear media. In this Section we briefly review the 
essential features of the theory of macroscopic quantum 
electrodynamics in nonlinear media, referring for a more 
exhaustive treatment to 0, 0, 0| . The standard ap- 
proach is the phenomenological quantization of the clas- 
sical macroscopic theory. The macroscopic description of 
the interaction of the electromagnetic field with matter 



is based on the polarization vector [17], [1^ [19 

Pi = x^^^ -E + x^^^ ■EE + x''^^ -EEE- 



X?,Ie,E,+x%E,E,Ei + ... {52) 

where E is the electric field, xtjli ('^ + 1) 

rank susceptibility tensor and the subscripts indicate 
the spatial and polarization components. For a lossless, 
nondispersive and uniform medium, the susceptibilities 
X^""^ are symmetric tensors. Considering only electric 
dipole interactions, the electric contribution to the elec- 
tromagnetic energy in the nonlinear medium is 



H = 



leoE^r,t)+Y,Xn{r) 



where 



Xn{r) 



EE...E 



n+\- 

In terms of the Fourier components of the field, 

. n+l 



(53) 



(54) 



n + l 



Y[dcj 



le 



1=1 

■x}''\u;2,...,u;n+i):EE...E . (55) 

In deriving Eq. (|55|l . it is tacitly assumed that x"('^) 
is independent on the wave-vector k (electric dipole ap- 
proximation) , that in the electromagnetic energy of the 
medium terms of order cudx^^^ /duj can be ignored, and 
that nonlinearities are small. 

The canonical quantization of the macroscopic field in 
a nonlinear medium is introduced by replacing the classi- 
cal field E{r, t) with the corresponding free-field Hilbert 
space operator 

-,1/2 



E(r,i) 



2epV 



i{k-r- 



(56) 



-n^ ft p-«(fc-''-'^r,A*) 



where V is the spatial volume, e^; ^ is the unit vec- 
tor encoding the wave-vector k and the polarization A, 
LUj^ — c\k\/n-i: is the angular frequency in the medium 
{n-^ X being the index of refraction), and a^^ is the cor- 
responding boson annihilation operator. Denoting by k 
the pair {k, A), and defining 



A, 



(57) 



the Fourier components of the quantum field are given 
by 



E(cj) = A]^{akeke''''M(tj - cjk) 



-ik-r 



6{UJ + LOU.)} 



(58) 
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The contribution of the n-th order nonHnearity to the 
quantum Hamiltonian can thus be obtained by replacing 
the Fourier components of the quantum field Eq. I|58() in 
Eq. 155|1 . Owing to the phase factors e"^*, many of the 
terms resulting from the Eq. (|55|l are rapidly oscillating 
and average to zero; they can then be neglected in the 
rotating wave approximation. The surviving terms cor- 
respond to sets of frequencies that, due to the constraint 
of energy conservation, satisfy the relation 

s n+1 

Xl'^ki = ^ t^ki , (59) 
and involve products of boson operators of the form 

and their hermitian conjugates. The occurrence of a par- 
ticular multiphoton process is selected by imposing the 
conservation of momentum. This is the so-called phase 
matching condition and, classically, corresponds to the 
synchronism of the phase velocity of the electric field and 
of the polarization waves. These conditions can be real- 
ized by exploiting the birefringent and dispersion prop- 
erties of anisotropic crystals. 

The relevant modes of the radiation involved in a non- 
linear parametric process can be determined by the con- 
dition Eq. H59f) and the corresponding phase-matching 
condition 

s n-fl 

J2k.= J2 (61) 

i—l i—s+1 

We will have to consider in detail the terms of the quan- 
tized Hamiltonian Eq. (|53|l up to the fifth-order contri- 
bution to the susceptibility tensor (n — 5), using the ex- 
pansion Eq. (|58|1 and assuming the matching conditions 
Eqs. and These terms are needed to generate 

the processes described by the Hamiltonian Eq. (|23|l and 
the associated HEMPSS, and are listed in Appendix IXI 

VIII. EXPERIMENTAL SETUPS 

A. Realizing the two— mode four— photon 
Hamiltonian 

We now specialize to a possible, but by no means 
unique, experimental setup the previous general results 
in the context of a multiple parametric approximation, 
in order to realize the particular multiphoton processes 
described by the two-mode Hamiltonians Eq. H23|l . 

We consider a nonlinear crystal, illuminated by twelve 
laser pumps at different frequencies : 

12 

=^S,e-*("-*-'''"-). (62) 

i=l 



We assume a lossless, nondispersive, uniform, noncen- 
trosymmetric medium, in order to exploit the properties 
of spatial and Kleinman symmetry of the susceptibility 
tensors, which can be taken to be real. Due to the high 
intensity of the lasers, all the pump fields are treated 
classically in this model. For simplicity, we will assume a 
single polarization of the light, and unidirectional prop- 
agating waves. Two quantum modes of the radiation 
field, respectively at incommensurate frequencies loi and 
UJ2, are excited in the parametric processes in the non- 
linear crystal. We impose the following relations on the 
frequencies of the laser pumps : 

Ql + fl2 = LUl + LU2 , n3+fli = 2(cJl + UJ2) , 

^5 + fie = 3^1 , rii + Qs = 30^2 , 

Qg + r^io = 2uJi — UJ2 , + r2i2 = 2ciJ2 ^ t^i (63) 

with the corresponding phase-matching conditions. With 
this choice, given lox and 0^2, each pair of frequencies 
is fixed by the corresponding phase matching condition. 
With a sufficiently large number of laser pumps, many 
related independent equations can be satisfied. Reduc- 
ing the number of pumps and of independent equations 
is certainly possible, at the price of leading to more com- 
plicated relations between the frequencies. 

The introduction of several pumps, with suitable fre- 
quencies, allows to study the multiphoton processes cor- 
responding to higher-order contributions to the suscepti- 
bility tensor. We will use conditions Eqs. 1)631) to obtain 
two-, three- and four- photon down conversion processes, 
and the associated mixed processes in the interacting 
Hamiltonian. We write all the contributions considered 
in our model in the interaction representation. Let us de- 
note by the coefficient of the generic interaction term 

a'^ 0^2 (^lO^ 1 describing the creation of j photons and the 
annihilation of I photons in the first mode, and the cre- 
ation of s photons and the annihilation of m photons in 
the second mode. The two-photon down conversion term 
due to third-order interaction reads: 

i?2 = K}^^a\a\E^E2 + h.c. , (64) 

with 

'^00 x^^H-^2;-t^i,-t^2,f^i)- 

The degenerate three-photon down conversion terms are 

Hsd = (KglJaf^s^e + h.c.) + (KO^af Sy^g + h.c.) , (65) 
with 

1^00 X^^H~^6; -(^1,-^1, -(^1,^5) , 

and 
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The semi-degenerate four-photon down conversion term 
is given by: 

H^sd = Kllal^al^EsE^ + h.c. , (66) 

with 

Kqo OC X*-^H-^4; -Wi, -Wi, -UJ2, -UJ2, ^3) ■ 

The semi-degenerate three-photon terms are 
Hssd = Klla\^a2EgEi() + h.c. , 
^3sd = Kllal^aiEiiEi2 + h.c. , (67) 

where 

«00 ^ x''^H-f^l0;-^^l,-Wl,W2,f^9), 

and 

i^m X^''H-^i2;-t^2,-W2,a;i,17ii) . 

In spite of their weaker contributions, we retain for con- 
venience the terms 

Hr — {kHo^^ a\aiEiE2 + h.c.) 

+ {Klla\al^a2EiE2 + h.c.) , (68) 

with 

nfl (X X^^H~^2; -t^l, -W2, -UJi,U!i,fli) , 

and 

'^Ol X^^\-^2; -(^1, -^2, -UJ2,UJ2, ^1) ■ 

FinaUy, we have to consider the Kerr contributions 

Hkerr = Kjooj^a^ + Ko2<^|^a2 -|- K];];a|a2aia2 , (69) 
where 

KqI OC X^^H-C^2; -W2, W2) , 

and 

In conclusion, the fuU interaction Hamiltonian reads 

Hj ~ H2 + H^d + -ff4sd + Hy,sd + -ffssd + Hr + Hkerr ■ (70) 

We stress again that our choice is only one of the many 
possible ones. One could envisage different configura- 
tions, using nonlinear crystals of particular anisotropy 
classes, and choosing different polarizations for the two 
quantum modes. 

Let us consider the interaction part of the Hamilto- 
nian Eq. (|23|l . We see that all its terms are included in 
Eq. H7UI) . and that the two Hamiltonians coincide if we 
impose the equality of the corresponding coefficients and 
assume 

'^20 ~ '^02 ~ . (71) 

Equality of the coefficients can be realized by tuning the 
complex amplitudes of the classical fields. 



B. Generating the four-photon HEMPSS 

Eqs. H44I) and 149|l suggest a possible method to re- 
alize experimentally the HEMPSS Eq. (|T7|l for the sim- 
plest, quadratic nonlinearity. Adopting the same criteria 
of Subsection IVIII Al we consider again a sample non- 
linear crystal, illuminated by eight classical pumps at 
different frequencies 

8 

^(+) = S,e-'(""*-'''--) . (72) 

i=l 

We impose the following relations 

r^i + il2 = 3ct;i , + il4 = 3^2 , 

il5 + fig = 2^1 — UJ2 , fir + ils = 2^2 — ^1 , (73) 

and coUinear phase- matching conditions. Wc then obtain 
an interaction Hamiltonian of the form 

Hj = KgOflf ^1^2 + ^lal'^EsEi 
+Kl°ia\^a2Er^E6 + n^^afaiErE^ + h.c. , (74) 

with 

'*oo x^^\~^2;-L0i;-i^i,-i^i,^i) , 



l^ol « X ' {-^i, ~t^2, -^2, -t^2, ^3) , 



«01 X^ '(-^6;-Wl,-Wl,W2,f^5), 



1^10 « X '(-f78;-t^2,-W2,t^l,f^7), 

where we have neglected the Kerr terms, due to their 
weaker contributions. 

Tuning the intensities of the external pumps, one can 
force the Hamiltonian Eq. (|74l) to coincide with the inter- 
action terms appearing in the exponent of Eq. (|49ll . The 
unitary operator Eq. (|49|l can thus be realized by the 
Hamiltonian unitary evolution e~^*^' , and the HEMPSS 
can be generated by imposing this unitary Hamiltonian 
evolution on a standard two-mode squeezed state. 

IX. CONCLUSIONS 

In this paper we define a canonical formalism for 
two-mode multiphoton processes. We achieve this goal 
exploiting generalized, nonlinear two-mode Bogoliubov 
transformations. These transformations are defined by 
introducing a largely arbitrary, nonlinear function of a 
heterodyne superposition of the fundamental mode op- 
erators, depending on two local oscillator angles (hetero- 
dyne mixing angles). The nonlinear transformations are 
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canonical once simple algebraic constraints, analogous to 
those holding for the linear Bogoliubov transformations, 
are imposed on the complex coefficients of the nonlin- 
ear mapping. The scheme generalizes the one introduced 
in the companion paper (Part I) jj . Among the possible 
choices of the nonlinear function we pay special attention 
to those associated to arbitrary n — th powers of the het- 
erodyne variables. This class of transformations defines 
canonical Hamiltonian models of 2n-photons nondegen- 
erate and degenerate processes. We have determined the 
common eigenstates of the transformed operators, and 
found their explicit form in the entangled state repre- 
sentation. These states form an overcomplete set and 
exhibit coherence and squeezing properties. They are 
non-Gaussian, entangled multiphoton states of bipar- 
tite systems. They can be defined as well by applying 
on the two-mode squeezed vacuum an unitary operator 
with nonlinear exponent that realizes a heterodyne mix- 
ing of the fundamental modes. We have thus named 
these states Heterodyne Multiphoton Squeezed States 
(HEMPSS). For quadratic nonhnearities, the HEMPSS 
realize a two-mode version of the cubic phase states . 
Thus the HEMPSS are potentially interesting candidates 
for schemes of quantum communication, and for the im- 
plementation of continuous variable quantum computing 
involving multiphoton processes. One of the most ap- 
pealing features of the HEMPSS concerns the study of 
their statistical properties. As in the single- mode case, 
but with even more striking effects, the statistics can be 
largely gauged by tuning the strength and the phases of 
the nonlinear interaction, and the local oscillator (hetero- 
dyne mixing) angles. For the case of lowest (quadratic) 
nonlinearity, we investigate possible routes to the exper- 
imental realization of canonical multiphoton Hamiltoni- 
ans, and the generation of the associated HEMPSS. Our 
proposal relies on the use of higher-order contributions 
to the susceptibility in nonlinear media, and on the engi- 
neering of suitable classical pumps and phase-matching 
conditions. 

In future work we intend to qualify and quantify the en- 
tanglement of the HEMPSS, and their possible use in the 
framework of the continuous variables universal quantum 
computation. The overcompleteness of the HEMPSS al- 
lows in principle to study the unitary evolutions ruled by 
any diagonalizable Hamiltonian associated to a quadratic 
nonlinearity, and by other Hamiltonian operators associ- 
ated to known group structures in terms of the trans- 
formed mode variables. In this way, we intend to intro- 
duce and characterize other classes of nonclassical multi- 
photon states amenable to analytic study. 



APPENDIX A 

a) X2 contribution 

We fix the frequencies , j "^ka such that = 



'^k2 + "^ka- Ignoring the oscillating terms, we obtain 

J^d\X2ir) - ^d^r J (_Q dc^ie-'"'*) 

•X^2)(cj2,t^3) : E(wi)E(a;2)E(w3) ^ 

2ieQ ^ AkiAk2Ak3X^^'(t^k2,t^k3) : eki^kaeks 

ki,k2,k3 

■aLak^flka / dVe-^C^i-'^^-''-) "- + /i.e. (Al) 
Jv 

Condition Eq. H61() eliminates the strong dependence on 
the phase mismatch Ak = ki — k2 — ka of the volume 
integral Eq. (|A1|) . 

The resulting nonlinear parametric processes (in a 
three wave interaction) are described by a Hamiltonian 
of the form 

H (xx^'^^a^bc + h.c. , (A2) 

where a, b, c are three different modes at frequency uja, 
LOb, u!c, respectively. Hamiltonian Eq. ljA2|l can describe: 
sum-frequency mixing for input b and c and LOa — Wb+Wc; 
non-degenerate parametric amplification for input a, and 
LUa = Ub+LOc; difference- frequency mixing for input a and 
c and uJb = uJc — (^a- 

If some of the modes in Hamiltonian Eq. (|A2|) degen- 
erate in the same mode (i.e. at the same frequency, wave 
vector and polarization), one obtains degenerate para- 
metric processes as : second harmonic generation for 
input b = c and coa = 2^;,; degenerate parametric am- 
plification for input a, and coa = 2ujb, with b — c\ other 
effects as optical rectification and Pockels effect involving 
d.c. fields. 

b) contribution 

We must now consider 

■^^\uj2,uj^,iOA) : E(wi)E(w2)E(lj3)E(u;4) . (A3) 

In this case, the relation Eq. H59|l splits in the two distinct 
conditions 

Wki = t^k2 + t^ka + t^k4 , (A4) 
^qi + ^q2 = ^qa + ^q4 > (^5) 

and the interaction Hamiltonian term takes the form 
H (X [Q;x'^H'^k2,^k3,^k4)aiLiak2ak3ak4 + h.c] + 

[^X^^^(-^q2>^q3>^q4)aqi"q2"q3aq4 + h-c] , (A6) 

where a, P are c-numbers, and Okj (i = 1, ..,4) are four 
different modes at frequency luu_^. 
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This term can give origin to a great variety of effects 
such as high order harmonic generation, Kerr effect etc. 

c) X4 contribution 

We manage now the term 



At last we consider the term 



1=1 

■^'^\uj2,uj'i,uji,uj5) ■■ E(tJi)E(cj2)E(cj3)E(cj4)E(cj5) . 

(A7) 

Also in this case the relation (|59|l gives two independent 
conditions: 



t^ki = t^ka + ^k3 + ^k4 + ^ks 



^qi + ^q2 = ^qa + ^q4 + ^qs > 



(A8) 



(A9) 



while the corresponding Hamiltonian contribution has a 
more complex structure with multimode interactions: 

H oc [(5x*-''''(t^k2,^k3,t^k4,^k5)aki°^k2ak3ak4ak5 



Jv 

^£0 J d?r j J|dw,e-*"^V^)(a;2,W3,W4,W5,W6) : 

E(c^i)E(cj2)E(w3)E(u;4)E(cj5)EK) ■ (All) 
In this case the relation lfCT|) gives three conditions: 

t^ki = t^k2 + ^ka + ^k4 + t^ks + t^ke , (A12) 



^qi + ^q2 = ^qa + ^q4 + ^qs + ^^qs 1 



(A13) 



^Pl + ^P2 + ^P3 = ^P4 + ^P5 + ^P6 ■ (A14) 

The corresponding Hamiltonian contribution is: 

H CX [ryix'^H^ka, Wkg, Wk4,t^k5,t^ke)aki°^k2ak3ak4ak5akB 
+V2X^^^ (-Wq2 , , Wq^ , Wq^ , Wq Jflj,^ flj^^ flqg aq^Oq^ Oq^ 



+7X^'H-'^q2,'^q3,'^q4,'^qJ««q3«q4«q.+/»-C.]. +,^3^(5) (-Cp,, -C.p3 , Cp,, C.p„ Ja^^ fl^^ at,3 flp, flp, flp. 



(AlO) 

Here (5, 7 are c-numbers, and Ok; (« = l,--5) are five 
different modes at frequency • 

d) contribution 



+ /I.C.] 



(A15) 



Here ry^ (i — 1, 2, 3) are c-numbers, and Okj (i = 1, 
are six different modes at frequency Wkj ■ 
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